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A  class  of  data  windows  is  proposed  for  uso  in  the*  estimation  of  the  power 
spectra  cf  stationary  stochastic  processes.  These  windows,  which  are  p.oneral i .'ra¬ 
tions  of  the  standard  Farzcu  filter,  me  constructed  by  the  use  of  appropriately 
normalized  B-spl  ir.es.  It  is  demonstrated  how  this  window  c  1  a.;s  nay  be  cc..  puta- 
tioraKy  impli  rented  using  c  Fast  Fourier  Transform  a  lunrf  thr.»,  I  he  efficiency  cf 
the  resulting  ptocedure  is  genera! lv  a  significant  improvement  over  the  stale  of 
the  art  vn  th  little  additional  computer  tine*  required. 


1.  INTROWT/frON' 


TM  5  psprr  f «s  concerned  with  the  problem  of  the  result  of  using  ilase  cfU  windows  is  a  riv,- 

digitally  estimating  the  power  spectrum  of  a  vide-  nlficar.t  inproveviunt  iu  estimates  of  power  spectra 
sense  stationary,  ert'odic,  stochastic  process  without  a  noticeehle  increase  in  the  commit  at » oral 

from  a  sample  function  of  finite  length.  Not  effort.  We  have  chosen  to  inc  continuous  notation 

long  after  Cooley  and  Tuhcy  {5}  introduced  the  throughout;  conversion  to  discrete  formulae  where 

fast  Fourier  transform  algorithm  (FFTl ,  its  appli-  appropriate  is  straight  forward, 
cation  to  the  estimation  of  power  spectra  was  dis-  2  PREI IN IN ARIFS 

cussed  by  Bingham,  Godfrey,  and  Tukey  f 2  j  and  by 

Welch  [18],  They  suggested  a  method  based  on  com-  We  first  outline  some  basic  notation  and  re- 

plex  demodulation  and  pointed  out  the  computational  suits  which  are  assumed.  Here  we  follow  Stein  ar;d 


these  cfiu  windows  is  a  sig¬ 


nificant  fnprovivient  ii'  (*st  i r.  si  r«  of  power  spoetta 
without  a  noticeehle  increase  in  the  commit  at i oral 
effort.  We  have  chosen  to  me  continuous  notation 
throughout;  conversion  to  discrete  formulae  whore 
appropriate  is  straightforward, 

2,  PRELIMINARIES 


We  first  outline  some  basic  notation  and  re¬ 


speed  advantages  which  can  be  obtained  by  using 
the  FFT  for  computing  Fourier  periodograms .  Welch 
proposed  a  direct  method  which  has  become  the 
standard  numerical  approach  to  the  problem  of 
powe**  spectral  estimation.  This  procedure  has 
three  parts:  subdividing  the  available  data,  using 
^  ^  a  data  window  t<’  compute  a  smoothed  spectral  csti- 
0  mate  ^or  *>ach  segment,  and  averaging  these  smoothed 

spectral  estimates.  The  present  paper  proposes 
{  ^  a  class  of  da fa  windows  derived  from  the  B -spline 
basis  functions  of  Schoenberg  (14J,  Empiric* llv, 


Weiss  (16).  I.1  denotes  the  Banach  space  of 

Lcbesgue  measurable  functions  which  arc  (Lcbesgue) 
integrable  over  the  real  line  F  .  For  f^L1  ,  the 
Fourier  transform  of  f  is  the  function  defined 
by 

f(c>  =  r„  f(t)C'2,Tiadt 

for  all  Z £ F  where  l  -  /-I  .  The  notation 
(...)*  denotes  the  Fourier  transform  of  (...). 

For  fjgCl.1  ,  the  convolution  h  -  f*g  Is  defined 
as  the  function  h (  hl  given  by 
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h(T>  «  ^  f(--t)g(t)dt 

for  .ill  ?6P,  Tf  ,  then  (f*g)'  «  fg. 

Similarly,  (fg)*  -  f**g  provided  the  functions  in¬ 
volved  are  such  that  all  operations  make  sense. 

For  a  >0,  fi4  represents  a  dilation  operator, 
i.e,,  \f(t)  *  f(at).  Whenever  f,g*k  L1,  it 

is  clear  that 

0.f*S.8>(->  =  l  *.(  (f*g)(?)l 

and  (S,f)*(?)  =  i«.v  f(t)  . 

We  shall  consider  a  sample  function  x(t) 

T  T 

on  a  finite  interval,  -  —  <  c  <  J  unless 
otherwise  stated,  of  a  stochastic  process 
(X(t)  : -«<  t  < «} .  We  assume  this  underlying 
stochastic  process  to  be  wide-3ense  stationary 
(Doob(8])  and  ergodic.  For  simplicity  we  take 
the  mean  value  to  be  zero.  Let  X(t)  have  the 
spectral  density  function  P(f).  The  terminology 
related  to  spectral  analysis  used  in  this  paper  is 
generally  that  of  Blackman  and  Titkey  *[3). 

It  is  apparent  that  high  resolution  and 
stability  (snail  bias  and  variance)  are  desired 
qualities  of  power  spectra  l  estimates.  However, 
a  compromise  voen  t:he  two  must  be  made,  fn 
order  to  i.btnin  a  stable  spectral  estimate  using 
the  classical  indirect  method,  an  appropriate 
function  of  lag  must  modify  the  sample  autcco- 
variance  function,  Fundanenta ily  the  indirect 
method  estimates  Pff)  via 

I’.  CO  -  f  C(t1t.(  r)e'2rlfT  dT 

i .  tn .  w  _  t 

wh.*re  l.(T)  is  the  finitely  supported  lag  window 
and  is  the  sample  aulocovar i ance  function 

t  ,(T-|t|)/2  _ 

i:(T)  f  — < r - 1  t  1  > /z  *(t+ jJx(e-T)<it. 

Usually  the  support  of  I.(r)  Ls  much  less  than 
the  entire  interval  -T<  r<  T  . 

The  arguments  for  using  lag  windows  (e.g., 
see  Jenkins  and  Watts  f 1013  ran  he  adapted  to 
show  that  a  procedure  which  utilizes  the  KFT  to 
estimate  power  spectra  directly  should  incorporate 
a  data  window.  That  is,  for  a  direct  procedure 
to  have  a  proper  balance  between  resolution  and 
stability,  the  data  itself  must  be  modified  by  a 
suitable  fut, cl  ion  of  time.  Thus ,  basically  the 


at  discrete  frequencies  where  W(t)  is  the 
finitcLv  supported  data  window. 

The  end  result  of  either  of  these  modifica¬ 
tions  is  that  we  in  fact  estimate  smoothed  values 
of  P(f).  More  specifically,  if  we  let  P„(f) 
represent  the  estimate  computed  by  either  the  in¬ 
direct  or  the  direct  approach,  we  have 

•*»VR[P#(f)l»  Q(f )*P(f ) 

where  the  average  may  be  taken  either  over  the 
ensemble  or  along  time  [3J.  <)(f)  is  called  the 

spectral  window  corresponding  to  the  window  used 
in  the  computation. 

The  direct  and  indirect  approaches  to  esti¬ 
mating  the  power  spectrum  of  a  stochastic  process 
are  not  fur, c’a, -.out a  1  ly  different.  Blackman  and 
Tukev  [31  have  derived  the  ccv.oct  i  c->  between  a 
data  window  Wft)  and  its  'h’e  m  valent"  lap,  win¬ 
dow  L(~) .  'bis  relationship  turns  out  to  he  the 
correlation  integral 

l.(-t  -  MflW(t,rMt  . 

By  "equivalent"  we  mean  that  th*_*  appropriate  use 
of  the  respective  "eqnival  -nt”  windows  yields 
power  spectral  estimates  which  ore  equal  in  ex¬ 
pectation,  Note  that  if  the  data  window  is  an 
even  function  of  time,  then  the  equivalent  lag 
window  is  proportional  to  the  convolution  of  the 
data  window  with  itself. 

At  this  point  it  should  be  clear  ll;.it  in  the 
indirect  method,  the  spectral  window  QO’)  cor¬ 
responding  to  the  lag  window  Lf--)  is  simplv  the 
Fourier  transform  of  Tf*"’!.  For  the  direct  method 
it  then  follows  that  the  spectra!  window  Q(f) 
corresponding  lo  the  data  window  W(t  )  is  given 
bv  Q(  f )  =  where  Iff-)  is  defined  to 

he  the  Fourier  transform  ot  'v(l).  J(f>  is 
called  the  frequenev  window  corresponding  to 
W(t).  We  again  refer  the  reader  to  Blackman  and 

fukey  (3). 
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3.  srur.r  data  » nwv/s 

1.1  MiVVI  VAT  ION 

lr  light  of  the-  rather  simple  rolat  tonsil  i  p 
between  a  data  window  apt!  its  iqu’v,ihiit  lav  win¬ 
dow,  one  should  not  ignore  the  »*>:tens  ivo  elicit 
which  has  boon  put  into  the  construction  of  la/ 
window*  hv  numerous  ant  hots  tior  the  properties  of 
Lhose  windows.  See  for  cyar.ple  the  cor??: 
given  hy  Jenkins  f  9  ?  or  by  larzen  [12,13). 

For  instance,  recall  the  continuity  c)a.?s  to 
which  these  functions  belong.  In  particular  we 
observe  that  one  of  the  lag  windows  which  has  been 
proposes!  by  Parzen  [13}  is  proportional  to  a  di¬ 
lated,  fourth  order  B-spline  basis  function.  One 
of  the  two  data  windows  suggested  by  Welch  [18] 
has  the. shape  of  a  second  order/  B^-spline  basis 
function;  he  recognized  that  his  proposed  window 
yields  the  Parzen  spectral  window*  Henceforth 
we  will  refer  to  this  window  which  is  given  by 


as  Inc  Parzen  data  window,  even  though  it  is  also 
a  member  of  the  family  of  spline  windows  proposed 
in  this  paper.  Today,  Welch's  procedure  (de¬ 
scribed  in  the  introduction)  incorporat ing  the 
Parzen  window  is  undoubtedly  the  most  widely  used 
method  of  direct  power  spectral  estimation. 
Bingham  et  al  (2}  proposed  the  following  data 
window  which  is  continuous  through  the  first  dc- 
rivat i ve: 


It  I <  .4T 


W(t)  =  < 


.At  <  It  I  <  .  vr 


^  0  ,  otherwise  . 

Wc  have  normalized  this  window  properly  although 
this  was  not  done  in  (2}. 

From  another  point  of  view,  we  note  that 
under  the  Fourier  transform  operator  the  effect 
of  an  occurrence  at  one  point  in  time  tends  to  bo 
spread  over  all  frequencies,  Furthermore,  time 
functions  with  corners  have  Fourier  transforms 
that  have  a  "ringing"  effect.  Thus  It  seems  de¬ 
sirable  to  have  a  data  window  of  a  high  continuity 


class  so  that  lh*‘  pm  I  i  f  i  c  ’  sample  junction  •Of*; 
snoot  lily  to  :'cr»i  at  it::  «  mis . 

Such  tot.s  i  c»»r  at  i  o~.s  r;ot  i  vat  or  e  t  o  exar  :r** 
the  p s  i  h  i  1  j  :  v  of  usi’v  data  witVi-ws  b«v»cd  upm, 
hi  gin- r  criui  l'i-.spli-.o  bps  if  i  e:.<  t  i ,  .  A  prolir.r- 

iiuii)  survey  is  cr.ro'.t  i  «v  i  «»*’,.  I'itst,  we  i  ay  cho%?s» 
the  dr  si  red  couti’.uity  class.  Stcondlv,  1 1. •  coi  r< 
ponding  spevt  ral  window*  have  th«  shape  ^  Hi — J 

so  that  the  magnitude  of  the  side  Jobe4*  rapidlv 
decreases  relative  to  the  rain  p»  .n:  h'.  the  contin¬ 
uity  (or  Older  k)  increases.  Tt  addition,  all 
spectral  windows  are  of  non -no  pat  i  ve  type  (Pfirzti* 

[  li  1  )  .  Thirdly,  deboot  [  7  ]  lias  mrived  a  genet  •;) 
n  1  porithm  which  can  be  used  for  Accurately  evalu¬ 
ating  any  order  li-spline  basis  function  with  nnv 
interval^cf  support  and  at.  any  nw  r.b  length. 

3.2  THE  ri’Nli/.MFdi lAh  Si'l  l NV  FFNCilONS 
Schoenberg  [  1'*  ]  introduced  the  spline  func¬ 
tions  in  1946  with  the  f cl  1  owing  def i ni t i on .  With 
k  a  positive  integer,  a  real  function  «A(l)  de¬ 
fined  for  I  (.  b  is  called  a  spline  fnreut'. 

of  order  k  or  degree  k-1  if  it  has  tk?  folli.w- 
ing  properties: 

(i)  A(t)  Is  of  class  C*"S(F). 

(ii)  «\(t)  is  composed  of  polynomial  arcs  of 
degree  at  most  k-1  . 

The  polynomial  arcs  arc  Joined  nowhere  other 
than  possibly  at  integers  n  if  k  is  even 
or  at  points  ti  +  J  if  k  is  odd.  The 
spline  is  said  to  have  knots  at  these  points 
If  we  define 

..k-l 


(lit) 


yr* 


c 


y>  0 


.  y<0 

l-i 

with  exception  of  the  special  case  0.  r  j  and 
let  6k  represent  Die  k"'  Older  central  ditti:- 
ence  operator  with  unity  step,  f.c.,  6~f(y)  * 

f(y)  and  6*f(y)  -  6l"!f(y+?  )  -  6‘~’  f(y-i  ) 
for  k  =  1,2,...,  then  the  function 

M*>  -  -^tV)7  ^t/1 

is  easily  shown  to  be  a  simple  c**.nmplr  oi  a  spline 

function  of  order  k.  M*(t)  is  called  the  funda- 
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mental  B -spline  basis  function  of  order  k  .  We 
list  a  few  of  the  many  important  properties  of 
Mk(t)  and  refer  the  reader  to  Schoenberg  [l'»,15]t 
Curry  anti  Schoenberg  f 6)  #  and  to  deBoor  [7J  for 

these  and  other  properties.  ’\(t)  is  positive 

k  k 

tn  the  open  interval  (-  r*,**  )  and.  ident ically 
^  1  k 

zero  elsewhere;  it  has  knots  at  -  j+v  where  v 
is  an  integer  satisfying  0  <  v  <  k.  Mk(c)  is 
an  even  function,  and  its  integral  over  R  is 
unity.  M*(t)  also  has  the  representation 


Mk(t)  * 


(k-1) ;  v=o 


k  k-V 

(-D 


•h» 


r- 


For  example,  th. 
by 

n 


e  fourth  order  function  is  given 

(2-!tl)3-f(i-|tl)3  ,  Itl  <  1 


M.(t)  = 


HIM1)3 


,t<  It  l<2 


,othorwi se. 
We  shall  require  two  integral  propertied: 


5.  co  - 


and  (Mp*M,)(t)  -  MpM(T>  . 

3.1  OAK  PROPOSED  CLASS  OF  DATA  WINTX.1V/S 
3.3.1  Ini reduction. 

The  class  of  spline  data  windows  we  propose 
consists  of  the  fundamental  spline  functions  each 
dilated  so  that  its  support  is  equal  to  the  length 
of  the  sample  function  and  normalized  so  that  the 
total  power  of  the  data  remains  invariant  when  the 
spline  data  windows  are  applied.  If  we  let  0* 
ieprese.it  the  normalization  factor,  then  wo  can 
use  the  results  of  .the  previous  sections  to  ob¬ 
tain  the  following  table  of  the  spline  data  win¬ 
dows  and  of  th;:  corresponding  frequency  windows, 
lag  windows,  and  spectral  windows. 

Since  the  total  power  >r  is  given  by 

-•  *  x*  (t)dt  *  T  Varfxj, 

‘T/r 


We  seek  constants  C*  depending  upon  the  order 
k  of  the  spline  data  window  such  that  also 


■  _v3 


where  we  use*  the  tilde  to  indicate  an  un -normal¬ 
ized  window.  The  oxprctai  ion  ot  this  equation  is 

<*  =  C;  i‘  '  ElU\(UlK[x3(01dC 
'-■•1 


=  C3  V .i r ( x  1  ’  W;(t)<lt. 

V-'/3 

Changing  the  variable  of  integration  we  obtain 
/ 

c;  =  k/f  M? (t)dt  . 

l2 

An  application  of  Parseval's  theorem  yields 

•*  •  -AW 


Therefore,  we  can  calculate  the  appropriate  nor¬ 
malizing  factor  for  any  order  spline  data  window: 
C}  “  1.0  and 


Ck 


2  V  L-llVlfc-O-'  c.t'1) 

i !  (2k-i - 1) 


for  kS?. 


Norma lining  fur  invariance  of  total  power  is 
equivalent  to  the  requirement 

i  (Oi  --  j,J!  Q(  f\df  -  i  . 


Data  Window 
W,(t) 

CdV«.(U  * 

Frequency 

W  i  ndow 

1,(0 

C,  f  (  .in(~ll7k)\k 
k  \~fl7k'  ) 

F.qu  i  valent 

Lag  Window 
U(t) 

7-[k- 6^,(0] 

i:3  /kO 

v-(T) 

Spectra l 

W  i  n<l  ow  ! 

%(i) 

1 

'”,1  /.iM("n/k)\  S‘ 

k  \  Tfi/k  ) 

vAiiij:  i. 
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We  point  out  that  Wx  is  simply  a  trunca¬ 
tion  function,  f.j  is  the  lay;  window  suggested 
hy  Bartlett  fl],  and  is  the  corresponding 

Hart  Lett  spectral  window.  Wa,  La,  and  are 
respectively  the  Parzen  data,  lag,  and  spectral 
windows , 

For  tho  problem  of  balancing  the  resolution 
and  stability  of  power  spectral  estimates,  ex¬ 
plicit  expressions  for  the  bias  and  variance  of 
smoothed  estimators  have  been  developed  in  the 
HteraJture.  Jenkins  and  Watts  [10]  and  others 
give  those  expressions  for  the  indirect  computing 
method.  Welch  1 17, 18]  derives  the  expected  value 
and  variance  of  the  discrete,  directly  computed 
estimates  for  a  stochastic  process  which  is 
ilaussian,  ia  addition  to  our  previous  assumptions 
about  the  process.  In  each  case  some  property  of 
i he  window  being  used  is* exp! ie i t ly  involved, 
five  practitioner,  however,  needs  seme  device  to 
aid  in  the  selection  of  a  window  that  will  yield 
?  compromise  suitable  tor  his  purposes.  Toe  most 
Cvuion  measure  of  windows  is  c tiled  bandwidth  (by 
j.vtlogy  with  filtering  problems)! 

Several  definitions  are  used  tor  bandwidth. 
Hackman  and  fukey  11)  use  ins  tod  the  term 
equivalent  width,  which  they  define  to  be  the 
ratio  of  the  square  of  the  Integral  of  the  spec¬ 
tral  window  to  the  integral  of  its  square.  If 
we  define 

i  P  tj®(n<n. 

J  -30 

t .hi n  the  equivalent  wUHi  is  given  by 

B,  (<))  }  (/!,,*}<  >). II 

tonkins  ami  Watts  110]  define  h  Midwidth  to  hr' 
the  width  of  the  rectangular  (spertral)  window 
fur  which  the  variance  ot  the  i.uU  recti  v  computed 
■'.p.*ttra1  estimate  is  the  same  .»•;  it  is  fw  the 
aiven  window.  Ibis  definition  of  bandwidth  \ie)d- 
i  M)| 

Bf<Q)  }  . 

At  has  been  noted,  tor  a  properly  normalized 
window  0(t)dt  1;  therefore  l>  ^  .  We 


mention  \\{  for  its  application  to  the  spline 
lag  windows  given  in  Inblc  1.  However,  its  de¬ 
fining  ptoju'it  .•  is  no  long*  i  tine  for  th-<  direct 
method  of  computation  which  we  advocate.  Fur 
the  class  of  spline  windows  we  have 


»«  “  .t. 


••l.’p  :  (-n'Ck  -  5  1  ) 

k'  ,  S,  ili-tk-i-l); 

thus  B ,  C Q ■„ )  is  easily  dotoraii  n,\i  by  invorting 

lk  .  By  using  least-squares  l.oclmiquus  Cor 
k  =  3, ,  }4,  wo  obtain  tho  upprnxi'nnt  ion 
II,  (QO  «  (l.bS(k)°  '1,:-  0.32) /T  . 

Parzen  (13)  doftnos  h.rvlw icltli  ,is  (In’  width  of  a 
root  any.  In  which  has  tho  srrao  aroa  ami  anno  r-.axi- 
ro.ii.n  hoigbt  as  tho  -j t  r  ’l  wi,.d,..\  1  !-.<-re - 

ru>',?  ’ '  nio.ir 

is.  lO,  — * _ _ 

f.iav  Q(f) 

and  Cor  Uu*  x;}l  in,  win.!,/,,".  wo  havo 

k3 

W  ,-.;7  • 

Another  commont  y  use!  :;V*as>ir  •  of  h  'idwidt  h  is*  the 
distance  \tuTui  the  ha  1  f -power  points  ot  tho 
spectral  window.  lh..;,  for  an  oven  aptciral  win¬ 
dow  cent -'fi  d  at  i  -  0,  H  ,(Q)  J  t  wlu-ro  7  is 
the  f  i  ►st  frequency  such  that  0(f)  1  <}  id).  T.et 
j  n  l  •  i  »•«  ,  res- :\t  rh,  spherical  Hos.„-l  f, met  ions 
of  t  he  t  i  rst  k  i  nd  .  I  hen 


Approx  i-  a  ions  f  n  •>.  ;:-ch>  :  :  r  s  :  -  •  1  with 

f’i'e.e  ■  ,  to  ;.iM,  ,  .....  ;  |  l  \\  .  I  ;  >  %  i 

quenllv  roll  usus  >  ue  :-ade  to  tin  ..-j  n  i --u 

botw-iu  ’h«  first  ;  .,i  ml  a.lovi. 

fhe  )•'■»  ol  (J,  ■  !  *  iv  ci!  it  t  t  wgi.vfv 

i 

n  Is  >  1  Uuj  f  I  4'  ’ :  b.-r . 

vat  i  >u  ;  .  ■  u  t  ho  i  >  h.i.i  ds*<.n  ?hi(  ike  1  i  as 
ot  'I'.etr  i*  f  i  siat  <  >rs  i;ur«  ius  as  the 

bandwidth  nl  li  e  wii.d.,v  ineres..s,  hut  that  the 
variance  Is  invrisid  •  p»  opor  I  i  ona  1  to  i  he  b.md- 
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.  dth.  U  can  also  be  argu  I  that  the  bnndv. jdth 
t  ;  1 1 .  window  rusi  be  at  least  as  small  as  the 
width  ‘-'l  the  iitMoVist  iwpor i •  at  detail  i:i  the 

»  t  I  Uvl. 

1  aMo  2  is  a  compilation  of  the  v  a  lions  con¬ 
st. tats  associated  with  the  spline  windows  given 
iii  i ;>:>*•'  1  with  order  up  t  o  k-21  . 


Order 

foil:'  t  ent 

Band 

;i  dt  l.ii  >  i 

k 

■ck 

(Q,> 

•  >•',  ‘Q  .  • 

1 

1 .PV:0 

l.soilo 

1. 

? 

1 .  7  0/0 

1.8393 

1 .3313 

1.2737 

2.‘.or.‘i 

2.301/ 

I .  f f  dd 

1.37!" 

A 

2  .  !•  ?  8  / 

2.  (.8  3/ 

1.91 V 3 

1. 8 >01 

3 

9  .SO', 9 

3.02(n 

2  .  1  37  1 

?  .  (J  <  •'  s 

,  6 

3. 90.' 7 

3.3219 

2 .  '  ,3,. 

a  p  '!  •,  )• 

7 

9. 37'0 

3 . 39  79 

2.3976 

2.9  1  (.s 

8 

9.8  398 

3.8336 

7.7379 

2  .  Srf  >1 

9 

3.2783 

9.0939 

2.0071 

'/  A  3  \ 

|t> 

S .  7  1 0.- 

A .  32  Od 

3.0669 

7.89  id 

11 

(1.1912 

/,  .  S3t)3 

3.2138 

1 . 0  1 1  1 

1? 

6 .  S9  2  7 

9.7923 

3.31.3** 

J.l/oi 

19 

6,97  :»8 

9,9399 

3.3030 

3.3011 

19 

7.3913 

3.1292 

3.6367 

3.9269 

IS 

7.7298 

3.3121 

3.7637 

3.3973 

If. 

8.1118 

3.9890 

3.8909 

3 . 6(.9  1 

17 

8. 9880 

3.  (>607 

9.0113 

3.7779 

18 

8.8387 

3.8263 

9.1286 

3.8873 

19 

9.2793 

S.QS82 

9.2927 

3 .  O'1'.  3 

20 

9 .  33  31 

6 . 19  36 

9.3338 

A . 00 8 A 

21 

9.9919 

6.2990 

A. A 621 

— 

A.  2000 

TABIE  2  . 

3.4  CONiTTINC  WITH  SDUNT.  DATA  WINDOWS 
Except  foi  the  choice  of  data  windows,  the 
method  we  propose  for  computing  power  spectra  is 
that  of  Welch  (18).  Departing  from  our  conven¬ 
tion,  we  assume  that  we  now  have  a  sample  func¬ 
tion  x(t)  tw'*r  the  interval  0<t^S.  Using 
for  example  a  bandwidth  criterion,  choose  the 
order  k  of  the  spline  window  to  he  used  and 
then  the  wind'*’  length  T.  fwMivfdc  y(t)  into 


N  segments  each  of  length  T  and  with  a  dis- 
pl:.r«iafiu  D  1-itu'icn  the  initial  points  of  teu- 
rmutl  ivc  seg::«  nl  s  .  Denote  the  S\  y  Milts  by  >;  «  ( t  )  . 


Thus  (N  *  I)!;  i  V  -  S,  and 

X;{t)  k|h  (j-  !)!>-(  lor  j  1 . N ;  -  t ;  '  . 

Welch  itst-s  the  chi'ice  of  I)  unlv  as  a  parameter 
to  increase  the  stAhiliry  of  the  resulting  esti¬ 
mated  power  spec n am.  In  the  c.in  ol  split,.  w  i  n  - 
do.;*!  an  additional  consideration  should  he  r.ade. 
11  I)  -  T ,  so’.ue  oi  the  data  is  weighted  substan¬ 
tially  heavier  than  the  remainder  .  As  can  be 
Seen  from  Figure  1,  the  inequity  is  mere  pro¬ 
ne.. need  as  the  order  of  the  spline  window  in- 
cren.es  in  the  sense  that  the  spline  data  win¬ 
dows  hoc ouie  prep.re*  si vel y  wore  peaked.  Since 
oik-  iu»rnally  assc-.es  that  all  portions  of  his 
data  are  equally  represent  at ive  of  the  stochastic 


process,  intuitively  the  order  of  the  spline 
window  should  also  influence  I>,  In  practice  a 
good  rule  of  thumb  for  relatively  short  s.v.plc 
functions  seems  to  he  p  =■  ••  T  for  k  -  A.  D  ^  *i' 


1  c' 
foi  k  -  10,  and  1)  ^  -,T  for  k-  ,":0  with  the 

ovei  lap  for  other  order  window:'  t.w.td  on  these. 

c 

Next  compute  the  smoothed  spectral  estimates 

P.co  *  =lf  x,(t)W.  (Oo*2mftdt  p,  j  -1 . N. 

1  ‘  -i/5  • 


Finally,  the  estimated  power  spectrum  is  given  by 

mo  ■  l  >■  mo  . 

w  i  - 1 


Regardless  of  the  order  of  spline  window 
used,  P,(f)  is  obviously  a  consistent  estima¬ 
tor  of  P(f)  (at  least)  whenever  the  stochastic 
process  is  Gauss  inn.  'lo  see  this,  the  express¬ 
ions  for  El P, (1)1  and  Var[rt(f)J  as  given  by 
Welch  [17,18]  need  only  be  compared  with  the 
convergence  properties  of  Lhe  spline  windows  as 
discussed  in  this  paper.  That  i ,  as  S  s*  vo 
have  the  following  type  convergence  for  a  sta¬ 
tionary  Gaussian  process:  if  T— a  such  that 

g -0,  then  E(  <P<  f>  -  H.(f))a)-0  . 

In  implementing  this  procedure  on  a  digital 
computer,  one  could  derive  a  closed  form  for  the 
K-'-jlfne  basis  functions  like  that  given  pre- 
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cov.irinnci*  fum-iion. 


v ious; ly  for  However,  for  a  general  p»*o- 

grfln  and  one  that  would  undoubtedly  bo  numeri¬ 
cally  preferable,  the  algorithm  of  deP.oor  (7]  is 
roe o»«me tided  for  evaluating  the  spline  functions. 

Jt  course,  by  symmetry  only  half  the  window  need 
be  directly  evaluated.  Brigham  [A]  thoroughly 
discusses  the  use  of  the  l*V-r  for  comput  ing  the 
discrete  Fourier  transform  which  approximates 
the  Fourier  integral  we  used  in  writing  Pj(f). 

4 .  EX  PE  R  l  MENTA  T.  RKS  U  LTS 

The  effectiveness  of  the  spline  data  windows 
is  readily  apparent  from  experiments  with  machine 
generated  autoregressive  processes,  power  spec¬ 
tral  estimates  computed  using  various  order  spline 
data  windows  are  compared  with  those  which  used 
the  Bingham -God f rey- Tukoy  and  Parzen  data  windows 
and  with  estimates  computed  by  the  Indirect  method 
:>f  dlnoknan  and  Tukoy  incorporating  the  Hanning 
lag  winder  [3J.  lhe  estimates  were  not  realized  bv 
dividing  out  the  virl.-rco  of  the  toel.jiie  pro¬ 
cess  ;  as  is  customary  when  one-sid«.d  power  spectra 
a**e  used,  the  magnitude  has  been  doubled. 

For  these  experiments  we  used  sample  func¬ 
tions  of  •uitorogress ive  processes  which  wore 
mac!)! no  generated  according  to 

3 

xi  Ojv‘.  f  «t  ,  t  l,...,lu24 

to 

where  ::i  is  the  order  of  the  process,  f«>t )  are 
the  regression  coefficients,  .and  e,  is  a  white 
noise  ptoCess  with  E{-\  |  U  and  Varls,)  -  1. 
i  in*  parti«;ilnr  process**;  which  wire  vised  in  tho 
illustrations  had  the  following  coot C jcients: 

(i)  First  order,  *j  -0.40 

(ii)  Second  order,  ».  1.00,  s ^  -0.50 

(iii)  Fourth  order,  a,  0.?j,  0.50, 

r,  -  -0.33,  ...  *  -0.50  . 

lhe  following  seller-,.'  i-h-nt  i  I  if*s  the  windows  in- 
vu!  /.-d  in  the  computation  of  the  various  power 
spec  t  ra  1  *•  .  t  i mates  . 

(i)  dinning  .\x;  estimate  computed  by  the  in¬ 
direct  thud  with  the  Hanning  lag  window 
{!)  and  using  xx  lags  of  the  sanplo  ;mto- 


(in  :w>-  XX;  Hin;g!..:.i -(.odfrey-Tuke/  ,I.J|  a 
window  was  used  vi  th  segments  of  tin* 
so,.:)],*  function  of  length  T  v;< 
samp les . 

U l > )  Parboil  ss:  1’aiv.on  ilat.i  window  was  used 
Willi  seppienf  s  of  tho  sa.:.plo  function  of 
lou:;Ih  T  -  xx  so.itplcs, 

(iv)  Spline  xx-yy-/./. ;  Spline  data  window  of 
wider  k  =  xx  w.T5  lined  wiih  se(>r-.ent  r  of 
Uio  i.v.-.ple  function  each  with  length 
T  J  yy  samples  end  displacement  D  -it 
samples  between  initial  points, 
ihe  graphs  i  1 1  ustrat irnj  these*  experiments  arc 
found  in  Appendix  1, 
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